Continuum states from time-dependent density functional theory 
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Linear response time-dependent density functional theory is used to study low-lying electronic 
continuum states of targets that can bind an extra electron. Exact formulas to extract scattering 
amplitudes from the susceptibility are derived in one dimension. A single-pole approximation for 
scattering phase shifts in three dimensions is shown to be more accurate than static exchange for 
singlet electron-He + scattering. 
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I. INTRODUCTION 

Ground-state density functional theory (DFT) 0, Q 
has become a popular electronic structure method in 
both quantum chemistry and solid-state physics, because 
modern approximations produce useful accuracy at mod- 
erate computational cost 0, Q . Now, electronic excita- 
tion energies of atoms and molecules are being calculated 
using linear response time-dependent density functional 
theory (TDDFT)0,0. In this scheme, bound — > bound 
transition energies are first approximated by the poles 
of the frequency-dependent Kohn-Sham (KS) density re- 
sponse function 0, 0, 0, and then corrected to the poles 
of the true response function, i.e., the true excitations. 
Bound — > continuum transitions, however, have not been 
treated in the same way because branch cuts of the KS 
and interacting response functions overlap, and because 
it is the phase shifts, rather than the energies, that are 
of interest in the scattering regime. 

Even though photoresponse was addressed in the early 
days of TDDFT 9] , and there has been a long interest of 
using density functional methods for the scattering prob- 
lem (e.g. ^3)> there is no formal theory based on TDDFT 
to study electron scattering. Such a theory might be par- 
ticularly useful in the emergent field of electron-impact 
chemistry in which large targets are struck by low 
energy electrons, so that bound- free correlations are sig- 
nificant [T^|. 

Several results relevant to this goal are presented here. 
First, we provide a proof of principle: the time-dependent 
response of an TV-electron ground state contains the scat- 
tering information for an electron scattering from the 
(N— l)-electron target, and this is accessible via TDDFT. 
Second, we show how this leads to practical ways of calcu- 
lating scattering phase shifts or, in one dimension, trans- 
mission amplitudes. Finally, in the simplest case, singlet 
scattering from He + , we find that TDDFT yields better 
results than static exchange, demonstrating its higher ac- 
curacy at low computational cost. 

Although we are not presenting a complete theory of 
electron scattering within TDDFT, such a theory can 
be built upon the rigorous results presented here, and 



become a competitive alternative to existing techniques 
for calculating electron-molecule scattering cross sections 
(e.g. Since continuum states are the current- 

carrying states in molecular electronic devices, we also 
anticipate applications of our 1-d results in the field of 
electronic transport through molecular wires |l4j . 



II. EXTRACTING SCATTERING 
INFORMATION FROM THE SUSCEPTIBILITY 

A. Theory 

Our starting point is the Dyson-like response equation 
that relates the susceptibility ^(r, 1 "';^) of a system of 
N interacting electrons with that of its ground-state KS 
analog, x s (r, r';u/) [g. In operator form (* indicates spa- 
tial and spin convolution): 



X = Xs +Xs * /hxc*X 



(1) 



where / HX c is the Hartree-exchange-correlation kernel 
(we use atomic units throughout): 
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a functional of the TV-electron ground-state density p(r) . 
In Eq.||5J|, v xc (r, t) is the time-dependent exchange- 
correlation potential induced when a time-dependent per- 
turbation is applied to the A^-electron ground state. We 
write the spin-decomposed susceptibility in the Lehman 
representation: 



with 



E ^(r)^(rO 

^ u-n n + io+ v 1 



(3) 



N 



F ncr (r) = <*o|Pa(r)|*») ; A,(r) = ^(r-r^ (4) 



2 



where is the ground state of the iV-electron system, 
^ n its n th excited state, and p a (v) is the cr-spin den- 
sity operator. In Eq.©, 0« is the — * transition 
frequency. 

For the remainder of this section, we restrict the anal- 
ysis to one dimension. Consider large distances, where 
the TV-electron ground-state density is dominated by the 
decay of the highest occupied KS orbital |l5j ; the ground- 
state wavefunction behaves as |l6[ : 



*o -> tpo \x 2 ,...x N ) 
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where tpo' 1 is the ground-state wavefunction of the (TV— 
l)-electron system (the target), Sq the spin function of 
the ground state and p(x) the TV-electron ground-state 
density. Similarly, 

<f„ -> ip^ 1 (x 2 , ...xn) ^j=- S n {a, a 2 , ...<jn) (6) 



where tp^ is an eigenstate of the target (labeled by 
n t ), S n is the spin function of the n th excited state, and 
<j)k n (x) a one-electron orbital. 

We focus on elastic scattering, so the contribution to 
Fncr (x) from channels where the target is excited vanishes 
as x — > oo due to orthogonality. Inserting Eqs.© and 
© into the Id- version of Eq. © , and taking into account 
the antisymmetry of both Wo and 



(x) — * \f p{x)(j)k n (x)5o,n t 
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The susceptibility at large distances is then obtained by 
inserting Eq.Q into Eq.©: 



X(X,X';UJ) = ^2xaa'(x,x' 
<j<j' 
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So,n t Ss ,s n + cc(u -> -oj) (8) 



"[R],[L]" implies that the integral is over both orbitals 
satisfying right and left boundary conditions: 



[R] 

CO*) 



e ±zfcx _|_ rke fikx ^ x _^ -pgQ 

t k e ±lkx , a: — > ±oo 



(10) 



When a; — > — oo and a;' = — x the integral of Eq.© 
is dominated by a term that oscillates in space with 
wavenumber 2a/2(e — 7) and amplitude given by the 
transmission amplitude for spin-conserving collisions t k 
at that wavenumber. Denoting this by x° sc , we obtain: 
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While this formula also applies to the KS system, its 
transmission t s (e) can be easily obtained by solving a 
potential scattering problem (i.e., scattering off the TV- 
electron ground-state KS potential). The exact ampli- 
tudes t(e) of the many-body problem are formally related 
to the t s (e) through Eas. (|ll|l and JJJ. This is the main 
result of this work: the time-dependent response of the 
TV-electron ground-state contains the scattering informa- 
tion, and is accessible via TDDFT. A potential scattering 
problem is solved first for the TV-electron ground-state KS 
potential, and the scattering amplitudes thus obtained 
are further corrected by / HX c to account for, e.g., polar- 
ization effects. 

While Eq. seems impractical as a basis for compu- 
tations, it leads to practical approximations. For exam- 
ple, if Eq.© is iterated once, we find through Eq.JTJJ the 
following useful distorted- wave-Born-type approximation 
for the transmission amplitude: 

t(e) = f s (e) + -^=((HOMo,e|/ H xc(£ + /)|HOMo,e)) , (12) 



2s 

where |homo, e)) is the product of the highest occupied 
KS orbital and the continuum KS orbital of energy s. 



Since only scattering states of the TV-electron optical po- 
tential contribute to the sum in Eq.© at large distances, 
it becomes an integral over wavenumbers k = \/2e, where 
s is the energy of the projectile electron: 
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In this notation, the functions 4>k n are box-normalized, 
and <f> kn (x) — (/>k(x)/y/L, where L — > oo is the length 
of the box. The transition frequency f2„ = E^ — Eq 
is now simply fl k = E^ 1 + k 2 /2 - E$ = k 2 /2 + I, 
where I is the first ionization potential of the TV-electron 
system, and Eq 1 and E^f are the ground and n th excited 
state energies of the M-electron system. The subscript 



B. Example 

We illustrate on a simple 1-d model of an electron scat- 
tering from a one-electron atom of nuclear charge Z [l7| 
in the weak interaction limit: 



i d 



i d 



2 dx( 



2 dx\ 



x 2 ) 



Electrons interact via a delta- function repulsion, scaled 
by A. With A = the ground state density is a simple 
exponential, analogous to hydrogenic atoms in 3d. 

(i) Exact solution in the weak interaction limit: First, 
we solve for the exact transmission amplitudes to first 
order in A using the static exchange method [3. The 
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FIG. 1: Real and imaginary parts of the KS transmission 
amplitude t s , and of the interacting singlet and triplet ampli- 
tudes (to first order in A), for the model system of Eg, 1131 , 
Z — 2 and A=0.5 in this plot. 



results for triplet (t tr i P ) and singlet (t S i ng ) scattering are: 

ik 



t 



trip 



to 



tn = 



Z + ik 



tsing — to + 2Xtl 



-ik 2 



(k-iZ) 2 (k + iZ) 



(14) 



(ii) Our TDDFT solution: The ground-state of the TV- 
electron system (TV = 2) is given to 0(A) by: 

1 S> (xio-i,x 2 o- 2 ) = -^=(j) (x 1 )(j)o(x2) [S ai fS a2 i - S ail S a2l ] , 

(15) 



where the orbital 4>o(x) satisfies [19J, 20]: 



1 d 2 

Y&x 2 ~ + X ^ ^ 

To first order in A, 

Mx) = VZe- z ^ + 



cf>o(x) = fMj)o(x) (16) 



2e- 3Z W + e- z W(4:Z\x\ -3) 

^ (17) 

The bare KS transmission amplitudes t s (e) character- 
ize the asymptotic behavior of the continuum states of 
v s (x) — —Z5(x) + \\(f> (x)\ 2 , and can be obtained 
to 0(A) by a distorted-wave Born approx. (see e.g. 
Ref.P): 



t Q + Aii 



(18) 



The result is plotted in Fig.l, along with the interacting 
singlet and triplet transmission amplitudes, Eas. l|14|l . 

We now apply Ea. HllI) to show that the / HX c-term of 
Eq.Q corrects the t s values to their exact singlet and 
triplet amplitudes. 

We need / HX c only to 0(A): 



where the / HXC of Eq.(JTJ is given to 0(A) by / HX = / H + 

/x = i /hx (= s/h hcrc )- E q-113 y ields: 

x{x,x'\uj) = Xs(x,x'\uj)+^- [ dx"xs{x,x";u)x(x" ,x';uj) 



(20) 

Since the ground state of the TV-electron system is a spin- 
singlet, the Kronecker delta 5s ,s n m Eq-© implies that 
only singlet scattering information may be extracted 
from x, whereas information about triplet scattering re- 
quires the magnetic susceptibility M. = ^2 aa i {po-'Yvaa 1 , 
related to the KS susceptibility by spin-TDDFT 22] : 



M(x,x';uj) = Xs(x,x';uj)- 



X 



dx"x s (x, x";uj)M(x", x'; u) 
(21) 



For either singlet or triplet case, since the correction 
to Xs is multiplied by A, the leading correction to t s (s) is 



determined by the same quantity, x^ *Xs ■> where Xs 1 ' is 
the th order approximation to the KS susceptibility (i.e. 

with v s (x) — vi°\x) = —Z5{x)). Its oscillatory part at 

, — 

large distances [23J (multiplied by yj p(x)p(—x) /ik, see 
Eq.(nj) is equal to \t\. We then find through Eas. ljTTj) . 
P0)l. and (EU) that 



>(o) 



tsing — t s 



Xt\ , ttrip — t s — Xti 

in agreement with Eas. (|14|) . 



(22) 



III. SINGLE POLE APPROXIMATION IN THE 
CONTINUUM 

We have yet to prove an analog of Eq. (|11|) for Coulomb 
repulsion in three dimensions. But here we use quantum- 
defect theory [24| to deduce the result at zero energy. 
Consider the I — Rydberg series of bound states con- 
verging to the first ionization threshold / of the TV- 
electron system: 



E n -E = I-l/ [2(n- Mn ) 2 ] 



(23) 



where fi n is the quantum defect of the n th excited state. 
Let 



-1/ [2(n-^ n f 



(24) 



be the KS orbital energies of that series. The true tran- 
sition frequencies Lu n = E n — Eq, are related through 
TDDFT to the KS frequencies Within 
the single-pole approximation (SPA) [||: 

u n = u) s>n + 2((homo, n|/ HX c(^n)|HOMo, n)) (25) 

Numerical studies j2^ suggest that Ap n — fi n — is 
a small number when n — ► oo. Expanding lu u around 
A/j, n = 0, and using / = -e HOMO , we find: 



(19) 



■ w Sj „ - A/i„/(n - fi s ,n) 3 
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FIG. 2: s-phase shifts as a function of energy for electron 
scattering from He + . Dashed lines: the line labeled KS corre- 
sponds to the phase shifts from the exact KS potential of the 
He atom; the other dashed lines correspond to the TDDFT 
singlet and triplet phase shifts calculated in the present work 
according to Eq. (l28ll . Solid lines: accurate wavefunction 
calculations of electron- He + scattering from Ref. 28]. The 
solid line in the center is the average of singlet and triplet 
phase shifts. Dotted lines: Static exchange calculations, from 
Ref. 29]. The asterisks at zero energy correspond to extrapo- 
lating the bound — > bound results of Ref.| 



bution). We also plot the results of a recent highly accu- 
rate wavefunction calculation (solid), and of static- 
exchange calculations [2^ (dotted). The results show 
that phase shifts from the TV-electron ground-state KS 
potential, S s (e), are an excellent approximation to the 
average of the true singlet/triplet phase shifts for an elec- 
tron scattering from the (N— l)-electron target, just as in 
our one-dimensional model; they also show that TDDFT, 
with existing approximations, works very well to correct 
scattering from the KS potential to the true scattering 
phase shifts, at least at low energies. In fact, for the 
singlet phase shifts, TDDFT does better than the com- 
putationally more demanding static exchange method, 
and for the triplet case TDDFT does only slightly worse. 
Even though Ea. (|28|) is, strictly speaking, only applica- 
ble at zero energy (marked with asterisks in Fig. 2), it 
clearly provides a good description for finite (low) ener- 
gies. It is remarkable that the antiparallel spin kernel, 
which is completely local in space and time, and whose 
value at each point is given by the exchange-correlation 
energy density of a uniform electron gas (evaluated at the 
ground-state density at that point), yields phase shifts 
for e-He + scattering with less than 20% error. Since a 
signature of density- functional methods is that, with the 
same functional approximations, exchange-correlation ef- 
fects are often better accounted for in larger systems, the 
present approach holds promise as a practical method for 
studying large targets. 



We conclude that, within the SPA, 

A^„ = -2(n-^„) 3 ((HOMo,7i|/ HX c(^n)|HOMo,n)). (27) 

Letting n — » oo, Seaton's theorem (jr liirin^oo [i n = 
S(e 0+)) 24] implies: 

5{e) = <5 s (e)-2^«HOMo,£|/ H xc(£ + /)|HOMo,e)) (28) 

a relation for the phase-shifts S in terms of the KS phase- 
shifts S s applicable when e — > + . The factor (n — ^t s ,n) 3 
of Eq. l|27|l gets absorbed into the energy-normalization 
factor of the KS continuum states. 

We illustrate in Fig. 2 the remarkable accuracy of 
Eq.|(SSJ when applied to the case of electron scattering 
from He + . For this system, an essentially exact ground- 
state potential for the N — 2 electron system is known. 
This was found by inverting the KS equation using the 
ground-state density of an extremely accurate wavefunc- 
tion calculation of the He atom 26]. We calculated the 
low-energy KS s-phase shifts from this potential, 5 s (e) 
(dashed line in the center, Fig. 2), and then corrected 
these phase shifts according to Ea. H28|) employing the 
BPG approximation to / HX0 27] (= adiabatic local den- 
sity approximation for antiparallel contribution to / HXC 
and exchange-only approximation for the parallel contri- 



IV. CONCLUSION 

To summarize, we have shown how, in one-dimension, 
scattering amplitudes may be obtained from TDDFT, 
and deduced the results for three-dimensions near zero 
energy for Coulombic systems. The ultimate goal is to ac- 
curately treat bound-free correlation for low energy elec- 
tron scattering from polyatomic molecules, with a com- 
putational cost lower than that of static exchange. An 
obvious limitation of the present approach is that it can 
only be applied to targets than bind an extra electron, 
and there is much work yet to be done: general proof of 
principle in three dimensions, testing of the accuracy of 
approximate ground-state KS potentials, developing and 
testing approximate solutions to the TDDFT Dyson-like 
equation, extending the methodology to cases where the 
anion has a sharp resonance rather than a ground state, 
etc. 
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